1/1 


/AD-A137  238 
UNCLASSIFIED 


NEUTRAL  SHEET  CURRENT  INTERRUPTION  A  F 1  ELD- AL 1GNED 
CURRENT  GENERATION  BY..IU)  CALIFORNIA  UNIV  LOS  ANGELES 
INST  OF  GEOPHYSICS  AND  PLANETARY  P..  T  SATO  FT  AL . 

MAR  83  AFGL-TR  04-0019  F  19628- 82 - K -00  1 9  F/G  4/1 


NL 


ty  4 


_ Unclassified _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGE 


la.  REPORT  SECURITY  CLASSIFICATION 

Unclassified _ 

3a.  SECURITY  CLASSIFICATION  AUTHORITY  *"  ” 

3b.  OECLASSIFICATION/OOWNGRADING  SCHEDULE 


REPORT  DOCUMENTATION  PAGE 

~  lib.  RESTRICTIVE  MARKINGS 


3.  OISTR I BUTION/AVAI LABILITY  OF  REPORT 

Approved  for  public  release; 
distribution  unlimited 


14.  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 


[S.  MONITORING  ORGANIZATION  REPORT  NUMBERISI 


N/A _ 

ia  NAME  OF  PERFORMING  ORGANIZATION 

UCLA 


Bb.  OFFICE  SYMBOL 
Inapplicable) 


Be.  ADDRESS  (City.  Suit  and  ZIP  Cod* ) 

Institute  of  Geophysics  &  Planetary 
Physics,  University  of  California 
Los  Angeles ,  CA  90024 


IBb.  OFFICE  SYMBOL 
inapplicable) 


Be.  NAME  OF  FUNDING/SPONSORING 
ORGANIZATION 

AFGL 


Be.  ADDRESS  (City.  Stmt *  and  ZIP  Cod*) 

Hanscom  AFB ,  MA  01731 


1 1.  TITLE  (Inelud*  Security  Claeti/ication)  Unc laSSl 


AFGL-  TR-84-0019 _ 

7a  NAME  OF  MONITORING  ORGANIZATION 

AFGL 

7b.  ADDRESS  (City.  Stmt*  and  ZIP  Cod*) 

Hanscom  AFB,  MA  01731 
Monitor /William  J.  Burke/PHG 

B.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 


10.  SOURCE  OF  FUNDING  NOS. 


PROGRAM 
ELEMENT  NO. 


62101F 


PROJECT 

TASK 

NO. 

NO. 

7601 

17 

(1  Current 

Generatio 

*  fnitw.ri  iMus.i r?  b  n 


i3.  personal authorisi  Three  Dimensional  Driven  Reconnection 
Tetsuya  Sato*.  Takaya  Hayashi*.  Raymond  J. Walker 


13a  type  of  REPORT 
Reprint 


is.  supplementary  notation  Reprinted  from  Geophysical  Research  Letters,  Vol.10,  No.3,pp  221-224, 
♦Institute  for  Fusion  Theory, Hiroshima  Univ, Hiroshima  730, Japan  March  1983 


13b.  TIME  COVERED 
FROM  _  TO 


■  bj  -Mb 


15.  PAGE  COUNT 

4 


COSAT  I  CODES  IS.  SUBJECT  TERMS  (Continue  on  revert*  if  necettary  and  identify  by  block  number ) 

- —————— - key  words:  neutral  sheet,  reconnection  model 

— - - current  interruption,  field-aligned  current 

20  03 


IB.  ABSTRACT  (Continue  on  revert*  tf  neceteery  and  identify  by  block  number) 

We  have  simulated  externally  drive n  reconnection  by  solving  the  magnetohy¬ 
drodynamic  equations  in  a  three-dimensional,  tail-like  geometry.  As  recon¬ 
nection  proceeds, the  tail  current  is  interrupted  locally  and  field  aligned 
currents  are  generated.  The  field  aligned  current  flows  towards  the  iono- 
shere  on  the  morning  side  and  away  from  the  ionosphere  in  the  evening.  The 
field  aligned  currents  flow  in  a  narrow  band  at  the  outer  edge  of  the  plasma 
sheet.  Thus,  the  simulation  demonstrates  that  the  nightside  substorm  curren 
system  is  a  natural  consequence  of  the  driven  reconnection  model . 


DIBTRIBUTION/AVAILABILITY  OF  ABSTRACT 
UNCLASSIPIEO/UNLIMITEO  S  same  AS  RPT.  □  OTIC  USERS  □ 


33a  NAME  OP  RESPONSIBLE  INOIVIOUAL 

WILLIAM  J.  BURKE,  Chief  yfr. 
Space  Plasm 


I4i  *  r?  ^  ru  ; 


21.  ABSTRACT  SECURITY  CLASSIPICAT 


Unclassified 


33b  TELEPHONE  NUMBER 
(Inmtmd*  Aim  Code 


nclasslfied 


T 


VI  : 


GEOPHYSICAL  RESEARCH  LETTERS,  VOL.  10,  NO.  3,  PAGES  221-224,  MARCH  1983 

AFCL-TT.  •  C  4  -  0  0  1  9 

Neutral  Sheet  Current  Interruption  and  Field-Aligned  Current 
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Abstract.  We  have  simulated  externally  drived  reconnection  by 
solving  the  magnetohydrodynamic  equations  in  a  three-dimen¬ 
sional,  tail-like  geometry.  As  reconnection  proceeds,  the  tail  cur¬ 
rent  is  interrupted  locally  and  field  aligned  currents  are  gener¬ 
ated.  The  field  aligned  current  flows  towards  the  ionosphere  on 
the  morning  side  and  away  from  the  ionosphere  in  the  evening. 
The  field  aligned  currentsflow  in  a  narrow  band  at  the  outer  edge 
of  the  plasma  sheet.  Thus,  the  simulation  demonstrates  that  the 
nightside  substorm  current  system  is  a  natural  consequence  of 
the  driven  reconnection  model. 

Introduction 

Field  aligned  currents  in  the  magnetosphere  have  been  the 
subject  of  intense  study  during  the  past  decade.  This  is  because 
they  are  believed  to  provide  us  with  key  information  on  mecha¬ 
nisms  coupling  the  solar  wind,  the  magnetosphere  and  the  iono¬ 
sphere.  Two  sets  of  large  scale  persistent  field  aligned  currents 
have  been  observed  within  the  auroral  oval  (Zmuda  and 
Armstrong,  1974;  Iijima  and  Potemra,  1976).  The  ionospheric 
Pederson  currents  which  close  these  field  aligned  currents  con¬ 
sume  sufficient  energy  (~  1011  watts)  that  the  field  aligned  cur¬ 
rents  would  disappear  within  a  few  minutes  without  a  continuous 
energy  supply  from  the  solar  wind. 

The  high  latitude  pair  of  field  aligned  currents,  region  1  cur¬ 
rents,  are  thought  to  be  related  to  polar  cap  convection  (Sato  and 
Iijima,  1979).  The  lower  latitude  region  2  currents  are  more 
variable  than  the  region  1  currents  and  are  believed  to  be  a 
diversion  of  the  ring  current  into  the  ionosphere  (Wolf,  1975;  Sato 
and  Iijima,  1979)  associated  with  the  divergence  of  the  gradient 
and  curvature  drift  current. 

In  addition  to  these  persistent  field  aligned  currents,  several 
observational  studies  have  inferred  a  localized  transient  field 
aligned  current  system  near  midnight  during  substorms 
(Akasofu  and  Meng,  1969;  McPherron  et  al.,  1973;  Bostrom, 
1974).  This  current  flows  into  the  ionosphere  on  the  morning  side 
and  away  on  the  evening  side  and  is  believed  to  connect  to  the 
ionospheric  westward  electrojet  (see  the  recent  reviews  by 
Kamide  ( 1982)  and  Baurqjohan  ( 1982)  for  a  complete  list  of  refer¬ 
ences).  Atkinson  (1966)  suggested  that  such  a  current  By  stem 
would  result  from  localized  reconnection  in  the  magnetotail  and 
Sato  (1982)  has  developed  a  model  for  the  resulting  field  aligned 
currents  based  on  his  2D  simulation  of  forced  reconnection.  Sato's 
model  has  been  sketched  in  Figure  1.  He  found  that  as  reconnec¬ 
tion  develops,  the  cross  tail  current  tends  to  concentrate  in  two 
thin,  slow  shock  layers  (Sato,  1979),  and  argued  that  field  aligned 
currents  will  form  at  the  shocks.  However,  to  verify  this  argu¬ 
ment,  the  simulations  must  be  carried  out  in  three  dimensions.  It 
should  be  noted  that  Akasofu  and  Kan  (Kan  and  Akasofu,  1978; 
Akasofu,  1980)  on  the  other  hand,  have  argued  that  reconnection 
is  not  necessary  to  generate  the  substorm  current  systems.  They 
have  proposed  a  current  interruption  model  in  which  the  cross 
tail  current  is  diverted  into  the  ionosphere  without  reconnection . 

Recently,  Bim  and  Hones  (1981)  have  presented  a  three- 
dimensional  (3D)  simulation  of  tail  dynamics.  Starting  with 
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Birn’s  (1979)  self-consistent  3D  tail  models,  they  solved  the  mag¬ 
netohydrodynamic  (MHD)  equations  throughout  the  tail.  In  this 
calculation,  the  tearing  mode  instability  was  driven  by  a  sudden 
increase  in  resistivity.  After  the  resistivity  was  turned  on,  the 
plasma  sheet  began  to  thin  and  an  X-type  neutral  line  formed. 
Bim  and  Hones  also  observed  field  aligned  currents  in  their 
model  substorm.  The  field  aligned  currents  were  directed  tail- 
ward  on  the  dawnside  and  toward  the  Earth  on  the  duskside.  The 
polarity  of  these  currents  is  opposite  to  that  of  the  substorm 
current  system  described  above.  They  suggested  that  their  cur¬ 
rent  system  is  responsible  for  the  region  2  currents. 

In  this  paper  we  study  the  generation  of  field  aligned  currents 
in  the  magnetotail  by  using  a  high  resolution  3D  MHD  simula¬ 
tion  of  externally  driven  reconnection.  In  particular,  our  results 
confirm  the  field  aligned  current  model  of  Sato  ( 1982 )  and  con¬ 
trasts  that  of  Bim  and  Hones  (1981 1. 

Numerical  Method 

In  this  calculation  we  used  a  rectangular  coordinate  version 
(MAG1C3R)  of  a  3D  MHD  simulation  code  which  was  originally 
designed  to  study  the  spheromak  tilting  disruption  in  nuclear 
fusion  plasmas  (MAGIC3C)  (Sato  and  Hayashi,  1982).  In  order  to 
maximize  the  spatial  resolution  of  the  model,  we  Btarted  with  a 
very  simple  neutral  sheet  configuration.  We  adopted  a  Harris 
(1962)  type  magnetic  field  configuration.  The  magnetic  field  and 
plasma  configurations  we  used  are 

B(x,y,z>  =  (B0tanh  (z/L),  0,0) 
p(x,y,z)  =  p„ 

T(x,y,z)  =  T„sech2<z/L) 

where  B  is  the  magnetic  field,  p  is  the  mass  density,  T  is  the 
temperature,  B„,  p„,  T„  are  constants,  and  L  is  the  half  width  of  the 
neutral  sheet  in  the  z  (north-south)  direction. 

We  assumed  that  reconnection  is  triggered  by  a  non-uniform 
compression  of  the  plasma  sheet  somewhere  in  the  tail  accompa- 


Fig.  1.  A  three  dimensional  model  of  externally  driven  reconnec¬ 
tion  showing  the  generation  of  field  aligned  currents  (after  Sato, 
1982).  The  dashed  line  shows  the  position  of  the  slow  shock  in  the 
noon-midnight  meridian.  Field  aligned  currents  flow  in  the  slow 
shock  region  for  y?0. 
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nied  by  the  inflation  of  the  magnetic  flux  in  the  lobe  region  (see 
Sato  (1979)  for  a  detailed  description  of  the  driven  reconnection 
model).  To  model  this,  we  expanded  the  2D  driven  reconnection 
simulations  of  Sato  and  Hayashi  (Hayashi  and  Sato,  1978;  Sato 
and  Hayashi,  1979;  Sato,  1979),  to  3D.  In  order  to  compress  the 
plasma  sheet,  we  injected  a  magnetized  plasma  (Poynting  flux, 
kinetic  energy  flux  and  thermal  energy  flux)  through  the  north 
and  south  boundary  planes  (z  =  ±  L,)  of  the  simulation  box. 
Initially,  the  incoming  mass  flow  pattern  was  such  that 

pv,  (z  =  ±  L,)  =  *  -^(1  +  cos^-)  (1  +  cos-t^-i  (1) 
4  L*  Ly 

where  v,  is  the  z-component  of  the  velocity,  L,,  Ly,  L,  are  the 
dimensions  of  the  simulation  box,  and  A„  is  a  constant.  pv„  and 
pVy  were  initially  set  to  zero.  Later  they  were  adjusted  so  that  the 
mass  flow  vector  was  always  perpendicular  to  B  at  the  boundary. 
The  other  boundaries  (x  =  ±  L,,  y  =  ±  Ly)  were  assumed  to  be  free 
boundaries  through  which  plasmas  can  freely  enter  or  exit.  The 
inflow  pattern  in  (1)  enabled  us  to  reduce  the  physical  domain  of 
the  simulation  box.  We  did  this  by  assuming  symmetry  (or  anti¬ 
symmetry  when  appropriate)  about  the  equator  <z=0),  the 
noon-midnight  meridian  (y=0),  and  the  neutral  line  <x=0).  Thus 
we  were  able  to  reduce  the  simulation  domain  by  a  factor  of  eight 
and  thereby  to  increase  the  resolution  of  the  model. 

Recent  particle  simulations  have  demonstrated  that  anoma¬ 
lous  resistivity  can  be  generated  by  the  lower  hybrid  drift  insta¬ 
bility  (Winske  and  Liewer,  1978;  Tanaka  and  Sato,  1981a)  and 
that  it  is  more  strongly  generated  and  lasts  much  longer  in  the 
presence  of  an  external  force  that  compresses  the  plasma  sheet 
(Tanaka  and  Sato,  1981b).  The  resistivity  is  roughly  proportional 
to  the  square  of  the  diamagnetic  drift  velocity  (Huba  et  al.,  1978). 
Thus  we  have  adopted  the  resistivity  (tj)  model  (Sato,  1979) 

t,  =  a(VD-Vc>2  for  VD  >  Vc 

where  VD  is  the  diamagnetic  drift  velocity  and  oc  and  Ve  are 
constants.  VD  =  J/ne  where  J  is  the  neutral  sheet  current,  n  is  the 
plasma  density  and  e  is  the  electron  charge. 

The  simulation  system  was  a  rectangular  box  with  dimensions 
L*  =  3L,  L,  =  5L  and  L*  =  2L.  This  was  implemented  on  a  41  x  40  x 
51  point  grid. 

In  the  actual  calculations,  all  variables  were  normalized  to  the 
following  parameters:  L,  length;  VA  *  B„/(p^)0)1/2,  velocity;  BOI 
magnetic  field;  pOJ  mass  density;  BJ^L,  current;  BQ2/2^,  pres¬ 
sure;  B2/p„.  energy;  p„LVA,  resistivity;  and  B„VA,  electric  field. 


Magnetic  Field  Flow 
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Fig.  2.  Contours  of  constant  vector  potential  and  plasma  flow 
vectors  in  the  midnight  meridian  plane  (y=0)  at  t  *  12.8  rA  and 
14.8  rA. 


Current  Density  (z=.36L) 


-202 

X(L)  Earth  -*■ 

Fig.  3.  Current  density  vectors  in  a  plane  parallel  to  the  equator  (z 
=  .36L)  at  t  =  14.6  xA. 


In  the  computer  run  presented  here  we  set  A,  =  .2,  a  =  0-02  and  Vc 
=  3  while  B„,  T„,  p„  and  VD  were  all  normalized  to  1. 

Tail  current  interruption 

In  Figure  2,  magnetic  field  lines  and  flow  vectors  have  been 
plotted  in  the  central  meridian  of  the  tail  (y =0)  for  two  times  (T  = 
12.6  rA  and  14.6  rA).  Prior  to  T  =  14.6  rA,  the  acceleration  rate  (not 
shown)  increased  with  time  while  after  T  =  14.6  rA  it  saturated 
and  remained  nearly  constant.  In  the  noon-midnight  meridian, 
the  results  are  very  similar  to  those  from  the  2D  simulation.  In 
particular,  a  slow  shock  structure  appears  and  the  plasma  is 
accelerated  to  approximately  the  local  magnetoeonic  speed. 

The  direct  analogy  between  the  2D  and  3D  models  is  invalid  in 
planes  parallel  to  the  equator.  In  Figure  3,  the  cross  tail  current 
density  has  been  plotted  in  a  plane  parallel  to  the  equator  and 
0.36L  above  it.  The  current  is  almost  completely  interrupted  in 
the  central  region.  A  substantial  part  of  the  current  is  diverted 
around  the  center.  The  x-component  of  the  current  earthward  of 
the  reconnection  region  is  towards  the  Earth  on  the  dawn  side 
and  away  from  the  Earth  on  the  dusk  side.  Since  the  field  is 
primarily  in  the  x-direction,  this  indicates  the  presence  of  field 
aligned  currents.  This  can  be  seen  more  clearly  in  Figure  4  where 
magnetic  field  lines  have  been  superimposed  upon  current  den¬ 
sity  vectors  in  a  meridian  plane  on  the  dusk  side  1.15L  from  y=0 
There  are  two  sets  of  current  vectors  of  interest.  First  near  the 
outer  edge  of  the  reconnected  field  lines  current  density  vectors 
with  large  field  aligned  components  are  observed.  The  field 
aligned  current  density  (J ,  =  4  •  fi/B)  is  20%  to  nearly  100%  of  the 
cross  tail  current,  and  as  we  saw  in  Figure  3,  the  field  aligned 
current*  near  dusk  flow  away  from  the  Earth.  These  field  aligned 
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Fig.  4.  Current  denaity  vectors  in  a  plane  parallel  to  the  noon- 
midnight  meridian  on  the  dusk  side  of  the  tail  (y  =  1.15L)  at  t  = 
14.6  rA.  Magnetic  field  lines  have  been  superimposed  on  the 
current  density  vectors. 


currents  are  observed  at  all  y  values  except  y  =0.  The  amplitude 
(J,)is  largest  at  y  =  1.67L. 

The  second  set  of  interesting  current  vectors  (Figure  4)  is  just 
equatorward  of  the  Held  aligned  currents.  Here,  the  current  flows 
earthward  and  toward  the  equator.  However,  this  large  earth¬ 
ward  current  density  does  not  indicate  an  earthward  field 
aligned  current  since  the  field  has  changed  direction  in  this 
region.  In  contrast  to  the  field  farther  from  z=0,  the  field  in  this 
region  is  mainly  north-south.  Most  of  this  current  is  normal  to  the 
field.  This  means  that  the  cross  tail  current  which  was  initially  in 
the  y -direction  now  has  an  earthward  component.  In  Figure  6  the 
current  density  vectors  have  been  plotted  in  the  equatorial  plane 
(*=0).  The  earthward  flow  can  be  seen  in  current  vortices  which 
have  formed  between  the  reconnection  region  and  the  region  of 
reconnected  field  lines.  In  the  region  of  reconnected  field  lines  the 
current  is  now  from  dusk  to  dawn  rather  than  dawn  todusk.  The  J 
x  £  force  in  this  region  thus  opposes  the  flow  from  the  reconnec¬ 
tion  region.  There  is  a  small  field  aligned  component  in  the  near 
equatorial  region.  However,  the  current  direction  is  still  outward 
from  the  Earth.  The  field  aligned  current  density  is  largest  in  the 
region  of  reconnected  field  lines  farthest  from  the  equator. 

That  the  field  aligned  currents  reside  on  the  outer  edge  of  the 
model  plasma  sheet  is  emphasized  in  Figure  6  where  pressure 
contours  have  been  superimposed  on  the  current  density  vectors. 
To  the  accuracy  of  the  model,  the  field  aligned  currents  occur  in 
the  sharp  pressure  gradient  which  characterizes  the  outer  edge. 
This  is  the  region  of  the  slow  shock. 

Discussion. 

Our  simulation  results  show  that  the  night  side  subetorm  cur¬ 
rent  system  is  a  natural  consequence  of  driven  magnetic  recon¬ 
nection.  In  our  model,  the  tail  current  is  locally  interrupted  by 
three  dimensional  externally  driven  reconnection  and  a  field 
aligned  current  system  is  generated  with  current  towards  the 
ionosphere  on  the  morning  side  and  away  from  it  on  the  evening 
side.  This  current  system  presumably  closes  in  the  westward 
eleetnqet.  We  should  point  out  that  we  have  not  modeled  the 
ionosphere  in  this  calculation.  The  ionosphere  was  treated  as  a 
perfect  conductor 

The  current  system  suggested  by  our  calculation  is  consistent 
with  the  current  wedge  model  moat  frequently  used  to  interpret 
both  magnetic  field  and  STARE  radar  observations  obtained  near 
midnight  during  suhetorms  (see  the  reviews  by  K amide  (1982) 


and  Baunyohan  (1982)  for  a  complete  list  of  references).  The 
present  results  contrast  the  field  aligned  current  system  found  in 
the  simulation  of  Bim  and  Hones  (1981).  In  their  model,  the 
currents  are  opposite  to  ours  ( i  .e. ,  away  from  the  ionosphere  i  n  the 
morning  and  into  the  ionosphere  in  the  evening).  However,  our 
results  do  not  necessarily  contradict  those  of  Birn  and  Hones 
since  the  models  use  different  boundary  conditions,  different 
initial  plasma  and  field  configurations  and  different  resistivity 
models.  Indeed  it  is  possible  that  mechanisms  analogous  to  both 
the  resistive  tearing  model  modeled  by  Birn  and  Hones  and  our 
driven  reconnection  mode  are  operating  in  the  tail. 

Since  the  field  aligned  currents  reside  at  the  edge  of  the  recon¬ 
nected  field  lines  in  our  model,  in  the  tail  we  would  expect  the 
substorm  associated  field  aligned  current  to  reside  in  a  narrow 
band  on  the  outer  edge  of  the  expanding  plasma  sheet.  Recent 
observations  from  the  ISEE  spacecraft  by  Kelly  et  al.  (1981)  are 
consistent  with  this  picture.  By  using  B  observations  from  both 
the  ISEE  1  and  2  spacecraft,  they  inferred  both  the  current 
direction  and  the  magnitude  of  the  current  density.  During 
plasma  sheet  expansions  they  observed  intense  field  aligned  cur¬ 
rents  limited  to  a  small  region  (500-2500  km  in  the  north-south 
direction)  near  the  edge  of  the  plasma  sheet.  The  current  direc¬ 
tion,  also,  is  consistent  with  our  model. 

The  field  aligned  currents  in  the  magnetosphere  can  be  calcu¬ 
lated  by  solving 
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Fig.  6.  The  same  as  Figu  re  4  but  with  contours  of  constant  plasma 
pressure  superimposed  on  the  current  density  vectors. 


and  J0  =  cB  x  VP/B2 
•  „  o  dv 


where  Jp  and  J,  are  the  parallel  and  perpendicular  current  den¬ 
sities,  B  is  the  magnetic  field,  p  is  the  mass  density,  (1  is  the 
vorticity,  <Jln  is  the  inertia  current,  N  is  the  density,  P  is  the 
pressure  and  v  is  the  velocity  ( Hasegawa  and  Sato,  1979;  Sato  and 
Iijima,  1979).  The  first  term  in  (2)  represents  a  discharge  current 
due  to  convective  motion  and  may  be  important  for  region  1 
currents  (Sato,  1982).  The  second  term  represents  a  diversion  of 
magnetospheric  current  due  to  a  magnetic  inhomogeneity.  This 
term  has  been  used  to  explain  the  region  2  field  aligned  currents 
by  diversion  of  the  ring  current  (Sato  and  Iijima,  1979).  The  third 
current  source  originates  where  there  is  a  density  gradient  in  the 
direction  of  the  inertia  current  and  is  usually  negligibly  small 
(Sato,  1982). 


Sato  (1982)  has  argued  that  the  second  term  in  (2)  is  responsi¬ 
ble  for  the  substorm  field  aligned  current  system  (Figure  1).  In 
the  driven  reconnection  simulations,  the  cross  tail  current  con¬ 


centrates  in  two  slow  shock  layers  (see  Figure  1 1  of  Sato,  1979).  In 
the  plasma  sheet  this  means  J  i  should  be  largest  at  the  outer 
edge.  As  reconnection  develops  the  field  intensity  at  the  slow 
shocks  increases  due  to  the  pile  up  of  field  lines  which  are  driven 
towards  the  neutral  sheet.  In  our  3D  reconnection  simulation  this 
pile  up  is  greatest  at,  y=0,  the  center  of  the  reconnection  region. 


$0 

In  the  slow  shocks,  — —  6,  points  toward  the  center  in  both  the 

dawn  and  dusk  sectors.  Thus  J ,  •  VB  >  0  in  the  dawn  sector  and 
J  i  •  VB  <  0  in  the  dusk  sector  From  (2),  the  field  aligned  current 
is  outward  from  the  ionosphere  on  the  dusk  side  and  inward 
towards  the  Earth  near  dawn.  Thus  in  our  simulation  both  the 
direction  of  the  current  and  the  concentration  of  the  largest 
currents  in  the  slow  shock  region  (Figure  6)  are  consistent  with 
the  predictions  of  Sato’s  model. 
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